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Abstract 

We calculate the structural, electronic, and magnetic properties of MnO 
from first principles, using the full-potential linearized augmented planewave 
method, with both local-density and generalized-gradient approximations to 
exchange and correlation. We find the ground state to be of rhombohedrally 
distorted Bl structure with compression along the [111] direction, antiferro- 
magnetic with type-II ordering, and insulating, consistent with experiment. 
We show that the distortion can be understood in terms of a Heisenberg 
model with distance dependent nearest-neighbor and next-nearest-neighbor 
interactions determined from first principles. Finally, we show that magnetic 
ordering can induce significant charge anisotropy, and give predictions for 
electric field gradients in the observed rhombohedrally distorted structure. 

PACS: 75.30.Et, 75.10.Lp, 75.50.Ee, 76.60.Gv 
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I. INTRODUCTION 



Since the recognition of their unusual insulating and magnetic properties over a half 
century ago, the first-row transition-metal monoxides have been the subject of much exper- 
imental and theoretical interest. Because they are highly correlated Mott or charge-transfer 
insulators,0H first-principles theoretical approaches based on density-functional (DFT) in 
the local-density approximation (LDA)i'S have proven inadequate for understanding certain 
properties such as electrical conductivity and magnetic moments, and new or extended ap- 
proaches such as self-interaction corrected (SIC), Hi orbital polarization corrected,! on-site 
Coulomb corrected (LDA+U)Jli optimized effective potential (OEP), and model GW0 
methods have been developed in order to address these deficiencies. Nevertheless, DFT-LDA 
based approaches have proven successful in understanding a number of ground-state proper- 
ties, most notably in the cases of MnO and NiO which are correctly predicted to be insulat- 
ing and antiferromagnetic with type-II ordering (though the gaps are underestimated) .HUlI 
Moreover, the use of the generalized gradient approximation (GGA)@ has been shown to 
yield further improvements. Here, we focus on the larger- moment compound, MnO, from 
the standpoint of DFT in both LDA and GGA approximations. 

Above the Neel temperature, T/v = 118 K, MnO is a paramagnetic insulator with the 
rocksalt (Bl) structure. Below T N , it is a type-II antiferromagnetic (AFII, Fig. [I]) insulator 
with rhombohedrally distorted Bl structure. The rhombohedral distortion takes the form 
of a compression along the [111] direction, taking cubic angles from 90° to 9O.62°.0 Based 
upon augmented spherical wave (ASW) LDA calculations, Oguchi et al. EE showed that the 
electronic structure is very sensitive to the magnetic ordering, and that the LDA can predict 
MnO to be an insulator, but only when the magnetic ordering is AFII, as it is experimentally. 
Based again upon ASW-LDA calculations, Terakura et a/.t3 subsequently explained the sta- 
bility of the AFII ordering relative to nonmagnetic (NM), ferromagnetic (FM), and type-I 
antiferromagnetic (AFI) orderings in terms of the strong dda coupling through the oxygens 
of opposite-spin cation sublattices in the AFII ordering vs. same-spin sublattices in the AFI. 
Using a linear combination of atomic orbitals (LCAO) LDA approach, Belkhir and HugelEl 
examined the splitting of d subbands, concluding that pd hybridization in the AFII structure 
enhances the splitting relative to that predicted by classical ligand field theory. Using a lin- 
earized augmented planewave (LAPW) approach, Dufek et a/.t3 showed that the GGA can 
yield further improvements to LDA results, finding a minor improvement of the gap from 
1.0 eV to 1.4 eV and of the Mn spin magnetic moments from 3.72fiB to 4.15/ie, relative to 
experimental values of 3.6-3.8 eV and 4.58-4. 79/ie, respectively. Using an LCAO Hartree- 
Fock (HF) approach, Towler et a/.0 calculated structural, magnetic, elastic, and vibrational 
properties. They determined the correct AFII ground state and a rhombohedral distortion 
angle of 90.47°, very close to the experimental value. Mackrodt et alM subsequently ap- 
plied the same methodology to the calculation of phase transitions under pressure. Using 
a linear muffin tin orbital atomic sphere approximation (LMTO-ASA) approach, Cohen et 
a/.il calculated low- and high-pressure properties of MnO in the context of an investigation 
of magnetic moment collapse, and found that the GGA yielded an improved lattice con- 
stant, bulk modulus, and magnetic moment relative to LDA results. More recently, using a 
planewave pseudopotential (PWPP) GGA approach, supplemented by LDA+U calculations, 
Fang et a/£3 investigated low- and high-pressure properties. They determined the correct 
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AFII ground state, a rhombohedral distortion angle of ~ 92°, somewhat larger than the 
experimental value, magnetic moments, and a phase transition to a metallic NiAs structure 
at high pressure. Massidda et a/.@ have investigated magnetic-order-induced anisotropies in 
linear-response properties using a combination of ab initio and model calculations. In the 
course of their study, they determined a ground-state distortion angle of 90.66°, very close 
to the experimental value, and showed that the zone-center optic phonon frequencies and 
Born effective charge tensor exhibit significant magnetic-order-induced anisotropies, even 
when assuming the perfect rocksalt structure, contrary to the assumption of cubic symme- 
try commonly made in interpreting experimental data. 

The mechanism underlying the rhombohedral distortion below T/v was examined by 
Kanamoriil in the late 1950 's. Based upon paramagnetic susceptibility data, it was con- 
cluded that there are significant exchange interactions between nearest-neighbor (NN) 
cations and, based upon Curie- Weiss and Neel temperature data, that these interactions 
are predominantly direct in nature. It was further concluded, consistent with earlier sug- 
gestions of Greenwald and Smart ,0 that the rhombohedral distortion below T/v is due to 
these NN interactions. Assumingdistance dependent NN and next-nearest-neighbor (NNN) 
interactions, Rodbell and OwencH subsequently used a molecular field approach to derive 
an expression for the distortion in terms of these interactions which supported the above 
conclusions, having explicit dependence only upon the NN interactions. Lines and Jones!! 
subsequently deduced the same expression from a Heisenberg model containing only NN and 
NNN interactions, and used a Green's function approaches in the random phase approxi- 
mation to evaluate the exchange constants based upon susceptibility data. More recently, 
Oguchi et a/@0 calculated exchange constants based on ab initio ASW-LDA calculations, 
employing the Korringa-Kohn-Rostoker coherent potential approximation (KKR-CPA). In 
this work they showed that interactions beyond NNN were in fact negligible, justifying the 
assumption made almost universally in previous work. The constants obtained, however, 
were about three times larger than those determined from experimental data. Among the 
reasons suggested for the discrepancy were the assumption that total-energy differences were 
well represented by eigenvalue-sum differences, and spherical approximations in the ASW 
and KKR-CPA calculations. Subsequently, Solovyev and TerakuraS calculated exchange 
constants based on LMTO-ASA-LDA calculations in the context of a discussion of single- 
particle theoretical approaches to MnO. In this work, they employed an expression^ for the 
exchange constants based on infinitesimal displacements from the ground state configuration 
rather than finite rotations in a random medium as in Refs. [T3| and ini and employed spher- 
ical approximations as before. The constants obtained were, however, again significantly 
larger than those determined from experimental data. The main reason suggested for the 
discrepancy was the underestimation of intra-atomic exchange and charge-transfer energies 
in the LDA. 

Here, we calculate the structural, electronic, and magnetic properties of MnO by the full- 
potential linearized augmented planewave method, using both LDA and GGA approxima- 
tions to exchange and correlation, thus improving upon spherical, pseudopotential, and/or 
strictly LDA approximations in earlier work. We show that the observed rhombohedral 
distortion can be understood in terms of a Heisenberg model with distance dependent NN 
and NNN interactions determined from ab initio calculations, without the CPA or spherical 
approximations. And lastly, we consider the charge anisotropy induced by magnetic ordering 
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through ab initio calculations of electric field gradients. 



II. CRYSTAL, ELECTRONIC, AND MAGNETIC STRUCTURE 

General potential LAPwil calculations were carried out for NM, FM, AFI, and AFII 
magnetic orderings. Well converged basis sets were employed with a planewave cutoff of 
18.9 Ry, corresponding to an average of ~ 109 basis functions per atom (the exact number 
being k-dependent).o Local orbital extensions were used for Mn I = 1,2,3 and O I = 
1, 2 channels to include higher lying semi-core states and relax linearization errors. Core 
states were treated fully relativistically in a self-consistent atomic-like approximation, while 
valence states were treated scalar relativistically. Increasing the planewave cutoff to 24.7 
Ry, corresponding to an average of ~ 159 basis functions per atom, yielded a change in 
total energy of less than 2.5 mRy/atom in LDA tests. Highly converged special k point 
setsS were employed consisting of 512 points in the rhombohedral (NM, FM, and AFII) 
Brillouin zones and 486 points in the tetragonal (AFI) zone. Doubling the number of k 
points in each direction yielded a change in total energy of less than 0.15 mRy/atom in 
LDA tests. LDA and GGA calculations employed the exchange-correlation functionals of 
Hedin and Lundqvist,@ and Perdew and Wang,0 respectively. As shown in Fig. the LDA 
calculations predict the ground state to be antiferromagnetic with AFII ordering, consistent 
with experimentJii and previous HF calculations^ and density-functional calculations^!!^! 
employing spherical or pseudopotential approximations. The GGA results are qualitatively 
similar, predicting the correct AFII ground state, but with a lattice constant of 8.38 a.u., 
much improved from the LDA value of 8.16 a.u. relative to the experimental value of 8.38 
a.u.S 

Figures |3| and f|(c) show the Brillouin zone and LDA band structure and density of states 
corresponding to the ground-state AFII ordering. The GGA results are qualitatively similar, 
but with a larger exchange splitting and gap of ~ 0.29 Ry and ss 0.105 Ry, respectively, 
compared to LDA values of ~ 0.27 Ry and ~ 0.075 Ry, due to lower occupied and higher 
unoccupied d states relative to O p states. In both cases, the large exchange splitting and 
narrow e g band are sufficient to produce an insulator, consistent with experimented and 
previous calculation Significantly, however, both LDA and GGA calculations 

predict other configurations to be metallic (Figs. £|(a),(b)), with the implication that the 
spin-disordered paramagnetic state would also be metallic, contrary to experiment. These 
results are consistent with the characterization of MnO as a correlated Mott-Hubbard insu- 
lator. 



III. RHOMBOHEDRAL DISTORTION 

Figure |5| shows LDA results for the total energy vs. rhombohedral strain in the ground- 
state AFII ordering at the experimental volume. The GGA results are essentially the same. 
The rhombohedral strain is determined by a parameter g through the volume conserving 
strain tensor 
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As shown in the figure, the calculations predict the ground state to be slightly compressed 
(g < 0) along the [111] direction. A least squares fit of the results yielded minimizing strains 
of g = —0.01456 and g = —0.01465, corresponding to deviations of 1.68° and 1.69° from 
the 90° cubic angle, for LDA and GGA calculations, respectively. A comparison with ex- 
periment and previous calculations is shown in Table |. We find that both LDA and GGA 
results overestimate the experimental angle somewhat, consistent with recent PWPP-GGA 
calculations^ but inconsistent with recent LAPW-LDA calculations^ On the other hand, 
Hartree-Fock calculations^ which produce larger gaps, underestimate the distortion some- 
what. Within simple tight-binding pictures, interatomic exchange interactions are inversely 
related to the gap (~ t 2 /I or ~ t 2 /U, where t is a hopping parameter and I and U are 
on-site interactions related to the size of the gap). As discussed below, the distortion is 
driven by the variation of NN exchange with distance, so the tendency to overestimate the 
distortion may reflect the incomplete treatment of Hubbard correlations and resulting small 
gaps associated with LDA and GGA approximations. 

The internal energy of the crystal can be partitioned into elastic and exchange parts. The 
elastic part is symmetric with respect to small rhombohedral distortions of the ideal cubic 
structure. Thus, it is sufficient to consider the exchange part to understand the mechanism 
of the observed compression along [111]. Within the model, the size of the distortion is just 
due to the competition of the exchange part which favors distortion and the non-spin-order 
dependent elastic energy, given in the harmonic approximation by the shear elastic constant, 
C44 for the non-spin-ordered case. To set the energy scale we calculated C44 = 66 GPa using 
the LDA at the experimental cell volume with a ferromagnetic ordering. Experimentally, 
the two terms cannot be readily separated, but above the magnetic ordering temperature a 
value of C44 near 75 GPa is measured with very strong temperature dependence in the region 
of We now show that the exchange part can be understood in terms of a Heisenberg 

Hamiltonian with distance dependent NN and NNN interactions determined from ab initio 
calculations. Let 

H=^ / J 1 S i -S j + £ J 2 S t -S r (1) 

NN NNN 

Considering a single formula unit, the above Hamiltonian reduces to 

H FM = (dJ 1 + 3J 2 )S 2 , (2) 



H AF1 = (-2J 1 + 3J 2 )£/ 2 , (3) 

and 

tfAFII = _ 3 J 2 1^ (4) 

where S is the magnitude of the cation spin, for FM, AFI, and AFII orderings, respectively. 
We deduce the exchange constants from the above expressions and ab initio total energies for 
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each ordering at the experimental volume, with S = 5/2. Our results are shown in Table ||, 
along with previous ab and semi-empiricaflS0 results for comparison. While 

there are still significant differences, our ab initio results agree in sign and order of magnitude 
with the semi-empirical ones, and are generally closer to the semi-empirical ones than pre- 
vious ab initio results employing the CPA and/or spherical approximations, implying better 
consistency with thermal and spin-wave data. Still, a significant difference exists, perhaps 
reflecting the role of Hubbard correlations, which might reduce the effective values of J via 
larger band gaps (recall that in the density functional calculations the ferromagnetic case 
is metallic). In any case, we find both NN and NNN interactions to be antiferromagnetic, 
consistent with the Goodenough-Kanamori rulesES and experiment. We note that since 
these calculations involve energy differences on the order of a few mRy, large k point sets 
were required in order to attain sufficient convergence. The above results were calculated 
using 4096 special k points in the full zone. Reduction of this set to 1728 points yielded 
differences of less than 0.003 mRy and 0.001 mRy for J x and J 2 , respectively, in LDA tests. 
Convergence with respect to planewave cutoff was also checked. Increasing the cutoff from 
18.9 Ry to 24.7 Ry yielded a difference of less than 0.002 mRy for both J\ and J 2 in LDA 
tests. Calculated Mn moments varied by less than 2.2% over all magnetic orderings in both 
LDA and GGA tests. 

We deduce the distance dependence of the interactions from calculations on distorted 
structures. For a small compression along [111], adjacent-plane NNs move closer together 
while in-plane NNs move farther apart, so that J\ — > J[ for adjacent-plane NNs while 
Ji — > J" for in-plane NNs; while for a volume conserving strain NNN distances remain 
unchanged to first order.0 Considering again a single formula unit, the above Hamiltonian 
(H) reduces to 

H FM = (3J[ + ZJ'{ + 3J 2 )S 2 , (5) 



H AFl = (-J' 1 -J? + 3J 2 )S 2 , (6) 

and 

F AFII = (_ 3J / + 3J // _ 3J ^ S 2 (?) 

for FM, AFI, and AFII orderings, respectively. We deduce the exchange constants as before 
from the above expressions and ab initio total energies corresponding to the appropriate 
magnetic orderings, at the experimental volume and distortion. Our results are shown in 
Table [TI], along with semi-empirical values for comparison. As the table shows, our ab 
initio results predict the same distance dependence as the semi-empirical ones: whatever 
the nature of the NN interactions, whether direct or indirect, their strength decreases with 
distance on the relevant scale, as may be expected. 

The observed compression along [111] now follows straightforwardly from the above es- 
tablished distance dependence and antiferromagnetic nature of the interactions. For a small 
compression angle 5 > 0, 

J[ = Ji(l + c5), 
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J'( = Ji(l-c5), 

and, for a volume conserving strain, J 2 remains unchanged to first order; where J\ > (an- 
tiferromagnetic NN interactions) and c > (strength decreases with distance) as established 
above. Thus, for the AFII structure, 

^strained ^unstrained ~~ J l c0 u j 

and so the compression lowers energy. Physically, for a small compression, the antiparallel- 
spin NNs move closer together while the parallel-spin NNs move farther apart, while NNN 
distances remain unchanged. And so, since the NN interaction is antiferromagnetic and 
decreases with distance, the energy is lowered. Thus, while decisive in determining the mag- 
netic orderingjl! the strong NNN coupling through the oxygens has nothing (to first order) 
to do with the observed rhombohedral distortion. Rather, it is the weaker antiferromagnetic 
NN coupling which is decisive. 



IV. MAGNETIC-ORDER-INDUCED CHARGE ANISOTROPY 

The electric field gradient (EFG) tensor, the second derivative of the Coulomb potential 
at the atomic position,^ provides a sensitive probe of the electronic charge distribution; and 
due to this sensitivity, very accurate distributions are required in order to calculate EFGs 
reliably. Blaha et a/@@ were the first to show that the LAPW method could be used to 
calculate accurate EFGs with no adjustable parameters. In subsequent worked it was shown 
that, due to the extreme sensitivity of the EFGs to inaccuracies in the density, the local 
orbital extension of the method (LAPW+LO),i! as we employ here, is sometimes required 
in order to obtain reliable results. 

Massidda et a/S have investigated magnetic-order-induced anisotropies of two nonmag- 
netic (spin-independent) linear-response properties of MnO, finding significant anisotropies 
in the zone-center optic phonon frequencies and Born effective charge tensor despite the 
fact that the atomic positions had cubic symmetry, so that the only noncubic ingredient 
was in the spin channel. Since strong Hubbard correlations are thought to separate charge 
and spin degrees of freedom in general, measurements of such couplings in comparison with 
density-functional calculations can be particularly illuminating. Here, we consider magnetic- 
order-induced anisotropy in a static, nonmagnetic quantity: the charge density itself. Since 
the EFG (the largest-magnitude eigenvalue of the EFG tensor) depends sensitively on the 
charge density and vanishes for cubic symmetry, it provides a sensitive measure of such 
anisotropy. And using this measure, our calculations show that magnetic ordering can in- 
duce significant charge anisotropy in MnO — on the order of that associated with noncubic 
crystal structure in some cases. Table [TT1] shows our LDA and GGA results for MnO in vari- 
ous magnetic and structural conformations. Cubic AFI and AFII structures show significant 
(up to ~ 1.3 x 10 21 V/m 2 ) EFGs, indicating significant anisotropy due solely to the magnetic 
ordering, with an order of magnitude difference in the Mn EFGs in the AFI and AFII order- 
ings. Generally speaking, GGA calculations are thought to provide a better description than 
corresponding LDA results for the ground state properties of materials which contain light 
elements and are not strongly correlated. For example, GGA equilibrium lattice parameters 
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and structural coordinates are frequently in better agreement with experiment than LDA 
results for such materials. However, while one may conjecture that this should also be the 
case for EFG's, this is not well established yet because the error bars on the experimental 
numbers are often larger than the LDA-GGA differences due to the uncertainty in the nu- 
clear quadrupole moments (note that in most cases relative LDA-GGA differences are not 
as large as in the present case). In fact, comparisons of calculated density functional EFG's 
with experimental quadrupole splittings were even used to obtain a more accurate value 
for the quadrupole moment of the Fe-57 nucleus.0 In strongly correlated systems, Coulomb 
correlations differentiate the orbitals more than LDA calculations. In open shell materials, 
this generally would lead to larger EFG's. For example, a strong effect of this type was 
demonstrated for the plane Cu site in YBa2Cu307.e3 On the other hand, in closed shell 
systems (like MnO), the effect may be opposite due to the suppression of charge degrees of 
freedom relative to spin, i.e. a lowering of the EFG. The last row of Table [HI] corresponds 
to the AFII structure with the observed rhombohedral distortion, and so constitutes our 
prediction for the observed EFGs. These values should be amenable to experimental mea- 
surement and, as mentioned, comparison with density-functional predictions should prove 
quite useful. 



V. SUMMARY AND CONCLUSIONS 



We have calculated the structural, electronic, and magnetic properties of MnO from 
first principles, using the LAPW method, with both LDA and GGA approximations to 
exchange and correlation. In both the LDA and GGA we found the ground state to be of 
rhombohedrally distorted Bl structure with compression along the [111] direction, antiferro- 
magnetic with AFII ordering, and insulating, consistent with experiment and previous HF 
and density-functional calculations employing spherical, pseudopotential, and/or strictly 
LDA approximations. The GGA was shown to yield a much improved lattice constant. 
However, highly converged LDA and GGA results were shown to overestimate the observed 
rhombohedral distortion, consistent with one recent PWPP-GGA calculation but inconsis- 
tent with others, including recent LAPW-LDA calculations. We showed that the distortion 
can be understood in terms of a Heisenberg model with distance dependent NN and NNN 
interactions determined from first principles, obtaining exchange constants in better agree- 
ment with semi-empirical values than previous ab initio results employing the CPA and/or 
spherical approximations, but still somewhat larger in general. Whereas the NNN coupling 
plays the decisive role in the magnetic ordering, the weaker NN coupling was shown to 
play the decisive role in the structural distortion. Finally, we showed that magnetic order- 
ing can induce significant charge anisotropy, on the order of that associated with noncubic 
structure in some cases, and gave density-functional predictions for EFGs in the observed 
rhombohedrally distorted structure. 
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FIG. 1. (a) AFI ([001]) and (b) AFII ([111]) magnetic orderings. 
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MnO 




FIG. 2. LAPW-LDA results for energy vs. volume of MnO, predicting AFII ground-state mag- 
netic ordering. 
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TABLES 



TABLE I. LAPW LDA and GGA rhombohedral distortion angles for AFII MnO, and compar- 
ison to previous calculations and experiment. 



LAPW 


PWPP-GGA 


LCAO-HF 


LAPW-LDA 


EXPT 


LDA GGA 


(Ref. 23) 


(Ref. H) 


(Ref. 24) 


(Ref. |18|) 



1.68° 1.69° « 2.0° 0.47° 0.66° 0.62 



TABLE II. Exchange constants (in K) for MnO based on LAPW LDA and GGA calculations, 
and comparison to other ab initio and semi-empirical results. 



■h 
■h 

A 



LAPW 
LDA GGA 



9.8 
24.5 
12.3 

7.4 



18.8 
33.0 
21.3 
16.4 



LMTO-ASA-LDA 



(Refs. Pljj 



24.5 
43.6 



ASW-LDA 
(Ref. |(]) 
30 
30 



Semi-emp 
(Ref. H) 



10 
11 



Semi-emp 
(Ref. |H) 
8^9 
10.3 
10.0 
7.9 



Semi-emp 



(Ref. gjj) 
8^5 
9.6 
9.9 
7.5 



TABLE III. LAPW LDA and GGA EFGs (in 10 21 V/m 2 ) for various conformations of MnO. 



The last row corresponds to the observed rhombohedrally distorted ground state. The asymmetry 
parameters vanish in all cases. 







LDA 


GGA 


AFI: 


Mn 


1.31 


0.65 




O 


0.18 


0.13 


AFII: 


Mn 


-0.13 


-0.05 




O 


0.22 


0.27 


AFII (distorted): 


Mn 


0.10 


0.20 




O 


-0.12 


-0.08 
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